Effect of superconducting fluctuations on ultrasound in unconventional 
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We study the renormalization of sound attenuation and sound velocity by fluctuation Cooper 
pairs in layered superconductors. We consider the influence of s- and d-wave symmetry of the 
fluctuating order parameter, on both longitudinal and transverse phonon modes. We show that 
both unconventional order parameter symmetry and transverse sound polarization suppress the AL 
and MT terms, while the DOS contribution is the least affected. The combination of these effects 
can change the sign of the overall fluctuation corrections above Tc. We also compare the results 
obtained using the Ginzburg-Landau formalism with a microscopic derivation of the fluctuation 
corrections to the sound velocity in both s- and d-wave superconductors. These calculations are 
motivated by ongoing ultrasound measurements in organic superconductors. 

PACS numbers: 74.25.Ld, 74.25.Fy, 74.40.+k 



I. INTRODUCTION 

Ultrasound is a very powerful tool to study the elec- 
tronic properties of metals. In systems such as Sr2RuO_ 
Q heavy fermions, and organic superconductors 
the ultrasound technique was used to probe the symme- 
try of the order parameter. In layered organic super- 
conductors, ultrasound attenuation was one of the key 
methods to determine the phase diagram. Q Compounds 
from the k-(ET)2-X family can exist in different phases, 
such as normal metal, superconductor or Mott insula- 
tor. Most phase transition lines can be associated with 
anomalies of the sound velocity and sound attenuation 
as a function of temperature or pressure. [3 In particu- 
lar, preliminary experiments indicate that the N-S phase 
transition is accompanied by a pronounced fluctuation 
region. This can be explained as an effect of transient 
Cooper pairs above Tc, or, equivalently, fluctuations of 
the superconducting order parameter. It is known that 
in layered organic materials the fluctuation region is es- 
pecially broad due to strong electronic anisotropy, low 
charge carrier concentration and relatively high Tc. The 
manifestation of superconducting fluctuations in various 
physical phenomena in organic superconductors has been 
experimentally confirmed. 0,0] 

The sound attenuation and the sound velocity in met- 
als are determined at low enough temperatures by the 
interaction of phonons with the electronic system. In the 
Frohlich model, phonons couple to the electronic den- 
sity, hence the electron-phonon interaction is momentum- 
independent. This approach is not valid, however, in 
many systems of interest. In the presence of impurities, 
one has to consider the electron system in a reference 
frame moving together with the ion lattice, 
deed, in an impure metal, the elastic scattering of elec- 
trons as well as perfect screening at small phonon mo- 
mentum and relaxation to equilibrium occur in this os- 
cillating frame. In this formalism, the electron-phonon 



coupling appears through the stress tensor, rather than 
through the density operator. 

In order to include polarization effects on sound prop- 
agating in highly anisotropic media such as organic su- 
perconductors, one has to consider a tight-binding model 
instead of a continuous one. In the tight-binding model, 
the electron-phonon interaction comes from the modula- 
tion of the hopping parameter induced by the lattice de- 
formation when the sound wave propagates through the 
crystal. The electron-phonon vertex then depends on the 
orientation of the phonon momentum and the phonon 
polarization. In such a model, where the stretching of 
specific crystal bonds induced by the lattice deformation 
contributes to the electron-phonon vertex, one can ex- 
plain for example that anisotropy in the attenuation of 
different phonon modes can be as large as few orders of 
magnitude in Sr2Ru04. 

In the superconducting state, the symmetry of the or- 
der parameter can also affect the temperature depen- 
dence of the sound attenuation. Namely, it does depend 
on the presence of nodes and on their orientation with re- 
spect to the phonon momentum and polarization. Thus, 
sound attenuation experiments can probe the type of su- 
perconducting pairing. 0, 

One would naturally expect that the superconducting 
fluctuations 0,0, Q will interact with ultrasound in a 
similar way, thus revealing the order parameter symme- 
try. The magnitude of the ultrasound renormalization by 
fluctuations is not necessarily small. Previously, the fluc- 
tuation corrections to longitudinal phonon mode propa- 
gating perpendicular to the conduction layers in quasi- 
2D material were found assuming s-wave symmetry of 
the energy gap. The model that uses the quasi-2D 
open Fermi-surface and the electron-phonon vertex orig- 
inating from the modulation of the electron interlayer 
hopping integral in the presence of the ultrasound wave 
predicts rounding of the sound attenuation and sound 
velocity temperature dependencies near T — s- Tc. 
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Since there is a substantial experimental interest in 
this subject, we present estimations for the fluctuation 
corrections to the sound attenuation and sound velocity 
for longitudinal and transverse ultrasound propagating 
perpendicular to the conduction planes, for s- and d-wave 
symmetries of the order parameter. 

The paper is organized as follows: in Section ^ we 
present a simple estimation for the fluctuation correc- 
tions to the sound velocity that can be obtained from the 
Ginzburg-Landau formalism in order to provide a simple 
phenomenological basis for the microscopic results. In 
Section IIIII we describe the microscopic model for the 
quasi-2D superconductor and for the electron-phonon in- 
teraction at various sound polarizations. In Section IIVI 
we review the sound attenuation and sound velocity in 
the mean field approximation (without superconducting 
fiuctuations). In Sections Ivl - IVII we give the supercon- 
ducting fluctuation corrections for various pairing sym- 
metries and phonon polarizations. Finally, we present a 
discussion of our results and of their relevance to exper- 
iment. 



II. GINZBURG-LANDAU APPROACH TO THE 
FLUCTUATION SOUND VELOCITY 



In this section, we obtain the fluctuation corrections to 
the sound velocity from the Ginzburg-Landau free-energy 
functional. Before that, we shall give, for completeness, 
a summary of elasticit y th eory that can be found in more 
details elsewhere. [T^ [13 Even though it is less detailed 
than the microscopic approach that is given later, the 
Ginzburg-Landau point of view allows one to develop 
physical intuition. It can also provide an independent 
check of the microscopic results that are given in the rest 
of the paper. 



A. Sound velocity jump at T = Tc 

The difference between superconducting and normal 
state free energies at temperatures close to Tc without 
fluctuations is 



AF. 



NS 



Fn — Fs 



8tt 



where a is the Ginzburg-Landau coefficient, and 
1 / dui dui 



'■^ 2 V dxj ' dxi 



(1) 



(2) 



is the (3x3) strain tensor that is defined through the 
displacement u at a point x. In a continuous elastic 
medium, obeys Hooke's law: 



(3) 



where Uij is the stress tensor and Cijki the elastic modu- 
lus tensor. Symmetry considerations allow us to reduce 



the number of independent components in Eq. (PJ, thus 
lowering the rank of and Uij and converting them into 
(1x6) vectors e,„ and Um with the following rule that re- 
lates the indices {ij} and {m}: 



yy 



yz 



xy 



6. (4) 



Similarly, the elastic modulus tensor Cijki is converted 
into a symmetric (6 x 6) tensor Cmn according to ji —^ m, 
kl n, and the rule 10}. The tetragonal system sym- 
metry reduces the number of independent components of 
Cmn to only 6 and imposes additional relations between 
them: cn = C22, C44 = C55, C12 = C21, C13 = C23 = C31 = 
C32, and all other off-diagonal components are 0. 

The equation of motion for sound waves results in the 
relation for the sound velocity Vg'. 



(5) 



where p is the density, and c = CijkiQiqjekei is the con- 
traction of the elastic modulus tensor with the phonon 
momentum and polarization vectors. In particular, for 
longitudinal phonons propagating along z-axis, the c be- 
ing expressed in terms of the (6 x 6) tensor is c = C33. 
For transverse phonons with the momentum q| |z and po- 
larization e lying in the xy plane, c = C13 for any angle 
If between the vector e and x-axis. 
Then, with the notations 



A. 



deaden ' 



(6) 



the discontinuity in the elastic modulus tensor at T ~ Tc 
that is related to the second derivative of AF^s with 
respect to strain is : 



ACrri.T). — ^ - 



1 92 AF, 
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T = Tr 



a^r r 



(7) 



where Vq is the unstrained volume. We shall consider 
the thermodynamics of the longitudinal phonons only, 
because for transverse modes vanishes by symmetry 
hence there is no discontinuity in Cm„ ■ Using the relation 
for the sound velocity Eq. (jSj), one finally obtains the 
sound velocity jump at T = Tc- 



Ac a^T^Tn 



2vsP SttVspVq 



(8) 



where the summation in Tm^n is done over all indices m 
and n that correspond only to non-zero components of 
c for the particular q and e. Finally, the sound velocity 
jump can also be related to the specific heat jump at 

T = Tr 



ACns — 



(9) 



as follows: llF 



Avs _ ACns 1 f I dTc 
Vs Tc 2p \ Vs dem 



^'^l (10) 
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This relation is quite general. It holds for both isotropic 
and anisotropic materials in the mean field approxima- 
tion. 



B. Fluctuation corrections to the sound velocity at 
T > Tc in the 3D isotropic case 

We can use the Ginzburg-Landau formalism to obtain 
the corrections to the sound velocity at T > Tc. In the 
presence of superconducting fluctuations, the free energy 
acquires the term : 



fi 



-T^ln 

P<PO 



(11) 



where e= {T — T^) /T , rj^ = l/{AmaT), and po ^ ^l\fnc 
is a cut-off. Then the most singular contributions are: 



^fl 



demder. 
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T{e + ijcP^y 



(12) 



Here the structure of the second term resembles the 3D 
fluctuation propagator L{^1 — 0,p) in the Gaussian fluc- 
tuation theory [191] , while the first term of Eq. (|12|1 looks 
like the fluctuation propagator squared. 

In Eq. H12|) . the main contribution is from the first 
term. If one changes summation to integration and sets 
the upper limit to infinity, in the limit e ^ it reads: 



po 



T(e + ry,p2)2 (2^) 



2V , 
2 JT^^dp = - 



(13) 



This expression is singular at T and clearly resem- 
bles the standard fluctuation conductivity contribution 
of the Aslamazov-Larkin type in the 3D case The 
second term of Eq. \12\ being integrated over momenta 
p is not singular as e — > 0. 

Thus, considering for simplicity only diagonal terms of 
the strain tensor, the corrections to the sound velocity 
are given by: 



1 
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Note that the standard derivation of the fluctuation 
corrections to the electronic specific heat, that requires 
taking second derivative of the thermodynamic potential 
with respect to temperature, results in the expression pol | 



1 



167r7^, 
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(15) 



that is in agreement with Eqs. H1U|) and (|14|) . The fluctu- 
ation corrections to sound velocity and specific heat here 
are connected in the same way as the mean-field jumps. 
This is not the case in the Lawrence-Doniach model, as 
we will discuss in the next section. 



C. Fluctuation corrections in Lawrence-Doniach 
model 

In the case of layered materials, one should use the 
Lawrence-Doniach quasi-2D free energy rather than the 
isotropic 3D one. At a low electron density, the quasi- 
particle spectrum can be taken to be 



i{p\hPz) 



p 
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where p|| denotes the intralayer component of the mo- 
mentum, pz the interlayer component, c the interlayer 
distance, and the interlayer hopping integral. The 
Fermi surface then has the form of a corrugated cylin- 
der with the one-particle density of states given by = 
m/{2T:c). Then, the fluctuation term analogous to Eq. 
(fTT|l reads: 



fl 
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P<P0 oi ■ 



\1 — cos 



(P.c)]} 
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Now 1 is related to the square of the in-plane correla- 
tion length, and the parameter r characterizes the ma- 
terial anisotropy. The modulation of hopping integral 
(5i_L by the strain due to the lattice deformation by 
the sound wave contributes to the electron-phonon in- 
teraction, that should be taken into account in the sec- 
ond derivative of Eq. p7() with respect to strain. 
Without going into all the details, we assume for the 



moment simple proportionality relations: 



t^. and 



5t± cx {dt±/dui){dui/dxj)5xj, where {dtj_/dui) con- 
tributes to the appropriate electron-phonon coupling con- 
stant g, and the combination (dui/dxj) can be rewritten 
in terms of the strain tensor Eq. Q. Therefore, one has: 



d^F, 
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(18) 



{r,„ - gt^T[l - cos(p,c)]}{r„ - gtiT[l - cos(p,c)]} 

+ + §[1 - cos(p^c)]| 

""Z" dp^p\\dp\\ - g^r[l - cos(p^c)] 

-^/c (27i-)2 e-H77||p2 4- §[1 -cos(p^c)]' 

We separate the expression in the last line of Eq. H18|l 
into two parts. The first one, without cos(p^c) in the 
numerator, results in the following contribution to the 
sound velocity: 



i^o(A„i„ - g^T) 



2mv 



sPVW 



In 



(19) 



The other part, with cos{pzc) in the numerator, gives a 
contribution of the type: 



4mi'spr7| 



(20) 
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In the integral in the second Une of Eq. (|18() . one 
can keep the terms in the numerator that do not contain 
elastic derivatives Tm- This leads to a contribution to 
the sound velocity of the following type: 



2mv, 



PV\\ 



1 



2e 
r 



-1 



(21) 



Keeping the Tm^n term in the numerator of the inte- 
gral in the second line of Eq. H18|l . and neglecting the 
hopping terms, gives the contribution: 



4:mVspTcr]\ \ \/e(e + r) 



(22) 



Finally, the cross-term that is proportional to both r„ 
and gt±^ results in the contribution: 



guojTm + r„) (yi - V^TT) 

2mvspr]\\ i/r(e -I- r) 



(23) 



that is less singular than Eq. H22() at small e and r. 

Note that the relation H1(J|) between the mean-field 
jumps of the sound velocity and specific heat at T ~ Tc 
is satisfied in the Lawrence-Doniach model. This comes 
about from the fact that in this model the difference be- 
tween superconducting and normal state free energies is 
still given by Eq. with the parameters modified for 
the quasi-2D case. On the other hand, the relation Eq. 
(|10|l is not satisfied for the fiuctuation corrections in this 
model. This is because we can neglect any T dependence 
of the interlayer hopping, hence it does not contribute to 
the specific heat. We will obtain later in the microscopic 
calculations in Section the expressions that have the 
temperature dependence of Eqs. ((T^ - lf^ . 



III. MICROSCOPIC MODEL 

A. Quasiparticle energy spectrum 

We use the model of quasi-2D square lattice with the 
in-plane period a, and the interlayer distance c > a. The 
quasiparticle energy spectrum in the normal state in the 
lattice model is given by: 



^(p) = -2t\i[cos{pxa) + cos{pya)] ~ p, 
— 2t± cos{pzc). 



(24) 



Here t|| and tj_ are the intralayer and interlayer hop- 
ping integrals. In highly anisotropic materials we have 
t\\ 3> t_L, and the electrons arc moving preferentially in 
the conduction layers. (In organic materials such as k- 
(ET)2Cu(NCS)2, for example, we can estimate the ratio 



i||/i-L ~ 4000 |2l|). Such a high anisotropy can raise the 
issue of incoherent electron motion in the perpendicular, 
z-direction. |22, l23l[ In the case of weakly incoherent in- 
terlayer motion, [23j the intralayer electron momentum is 



conserved in the tunneling process and the electron wave 
function in adjacent layers has some overlap, but there 
are many in-plane collisions between tunneling events. 
In this well as for coherent interlayer motion, the 

use of the quasiparticle spectrum Eq. H24() is legitimate. 
In the case of strongly incoherent interlayer transport, 
however, the intralayer electron momentum is not con- 
served in the tunneling processes between adjacent lay- 
ers because the tunneling can be accompanied by strong 
elastic or inelastic processes. This case needs a separate 
treatment. 

The tight-binding spectrum Eq. H24|) for app 1 
reduces to the form Eq. (|16|l . In most cases considered 
here the use of the low-density spectrum Eq. H16|) in the 
microscopic calculations is well justified. 



B. Electron-phonon vertex 

The electron-phonon part of the Hamiltonian, accord- 
ing to Walker, Smith, and Samokhin 0, is given by 



e—ph 



-V2^ y Gn(^^ 



k,p,R 



1/2 



(25) 



(^k-R^ (^-Rj (cosp-R) cj,_^k,<TCp,CT (a^k + «k 



where k is a unit vector in the direction of the phonon 
momentum, e a unit vector for the phonon polarization, 
R are the nearest bonds that are stretched by the sound 
wave. We have also defined uiq (k) = Vgk the sound fre- 
quency, Vs the sound velocity, M the ion mass, N the 
number of unit cells, Gr a constant that depends on 
the derivative of the hopping integral along the bond 
R with respect to the strain, and finally ajj-' and Cp^t 
are, respectively, destruction and creation operators for 
phonons and for electrons of spin a. 

For longitudinal phonons propagating in the perpen- 
dicular direction, the main terms in the sum over R in 
Eq. H26|) are those that contain the nearest neighbor in- 
terlayer bonds R = ±c (see Fig. ^a)). For these bonds, 
the wave vector and the phonon polarization will satisfy 
(k-R)(e-R) = 1. The electron-phonon vertex becomes 



ripHp) = 9L cos(p^c). 



(26) 



where is a constant. 

For transverse phonons, stretching of the nearest 
neighbor bonds Fig. ^^a) does not contribute to the 
sound attenuation as for these bonds (e-R) = 0. In other 
words, the shear wave does not interact with electrons to 
this order. One should consider next-to-nearest neighbor 
bonds, that in our particular case are the diagonals of 
the sides of 3D a x a x c unit cell. These bonds belong to 
the set {(±a, 0, ±c), (0, ±a, ±c)}. Let us assume that the 
phonon polarization makes an angle if with the x-axis in 
the plane, such that e = {costp, sint/j, 0}. The summa- 
tion over the next-to-nearest neighbor bonds results in 
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±a±c 




FIG. 1: (Color online) Side view of the basis vectors that de- 
termine the electron-phonon vertex for (a) longitudinal (L), 
and (b) transverse (T) phonons propagating perpendicular to 
the layers. Here a is the in-plane lattice period, c the inter- 
layer distance, and R are the nearest neighbor distances for 
the bonds that give the leading contribution to the electron- 
phonon vertex, for both types of polarization. 



an electron-phonon vertex of the form 

rir(<(5, p) = gT[sm{pxa) cos(p + sm{pya) sinip] sin(p^c), 

(27) 

where gr involves the derivative of the transverse hop- 
ping tx with respect to strain. In the particular case 
of sound polarization parallel to the x-axis, the nearest 
neighbor bonds that give the leading contribution to the 
sound attenuation, are determined by: R = ±a ± c (see 
Fig. nib)). Summation over these bonds yields 



rip (</5 = 0, p) = gr sin(p^a) sin(p^c) 



(28) 



In a similar manner, for phonon polarization oriented 
along y-axis, the vertex is 

rip(^ = 7T/2, p) = gr sinipya) sin(p,c). (29) 

It is natural to expect that the transverse coupling gT 
is significantly less than the longitudinal one g^. Still, 
it can provide a sizable effect, according to preliminary 
experiments in organics that suggest that the mechanism 
of the sound attenuation has electronic origin at T ^ Tc. 
At the same time, we neglect the contributions from all 
further neighbor bonds that lie for example along the 
spatial diagonals of the 3D unit cell, assuming that the 
magnitudes of corresponding hopping integrals decay fast 
enough. 

Finally, note that the single-particle spectrum Eq. 124|) 
should be modified because of the diagonal hopping tx 
mentioned above. That will introduce the term given 
by Eq. 158() . However, we will see in the next section 
that this effect leads to small corrections, even in the 
Aslamazov-Larkin diagram, where it is potentially im- 
portant, because tx/tj_ <^ 1. 



IV. ULTRASOUND ATTENUATION AND 
SOUND VELOCITY 

A. Normal state 

We start with the sound attenuation and sound ve- 
locity in the normal state. The sound attenuation coef- 
ficient is determined by the imaginary part 7(fc) of the 
complex frequency uj{k) where the pole of the phonon 
Green function D(k, w) is located. This quantity obeys 
Dyson's equation [2J| 



-D-i(k,w,) = [Z?°(k,cj,)] '-^n(k,c^,) 



(30) 



Expressed in bosonic Matsubara frequencies u!^ — 2tivT 
using units feg = 1 7i = 1, with v an integer, the quantity 
i^"(k,^.) 



DO(k,( 



-^o'(k) 



(31) 



is the phonon propagator in the non-interacting case, and 
n(k,a;j/) is the phonon self-energy: 



(32) 



j:j^rlpG{en,p)G{en +uju,p + k). 



In Eq. 1)32(1 the electron Green function G(p, e„) at finite 
temperature in the presence of impurities, is given by 



G(p,£n) 



1 



e(p)' 



(33) 



where ^(p) is the quasiparticle energy and where we de- 
fined £„ = £„ -I- l/(2T)sign(£„) with £„ = ■nT{2n+ 1) the 
Matsubara frequency and r the electronic elastic scatter- 
ing time. The quantity z/(2T)sign(£„) is the imaginary 
part of the quasiparticle self-energy. The real part of the 
self-energy is constant and is absorbed in the definition 
of the chemical potential. 

The scattering time r in Eq. (|32|l is the only effect of 
impurity averaging in the polarization operator. We do 
not include the corrections to the electron-phonon ver- 
tices Fgp. In the case of the current-current correlator, 
vertex corrections lead to the replacement of the scatter- 
ing time T by its transport analog Ttr- |24j For s-wave 
scattering, vertex corrections vanish because the vector 
vertex ev^ averages to zero upon angular integration at 
vanishing external momentum k. For stress tensor cor- 
relator, that we need in tight-binding limit, one can use 
results of Schmid 0. He showed that in the continuum 
limit, taking into account perfect screening, there is no 
impurity diffusion enhancement of the electron-phonon 
vertex in the case of transverse phonons, and that for 
longitudinal phonons this effect is negligible in the hy- 
drodynamic limit fc^ <JC 1 and <x'r<Cl. Physically, this 
comes about from the fact that the calculation should be 
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done in the moving frame and that screening is perfect at 
long wavelengths. In our case, the analogous argument of 
electroneutrality leads to the replacement of the stress 
vertex F hy F — (F), where the average accounts for the 
chemical potential shift. [2^ That average is precisely 
what is needed to make the impurity vertex correction 
vanish. In addition, in our specific case, (F) = to the 
order in phonon wave vector k that we need. 

As usual, after integration over p and summation over 
e„ in Eq. (|32|l . one should make the analytic continua- 
tion of the external phonon frequency following the rule 
icoi, — > Lu + iS. For sound propagating along z-axis, in 
the hydrodynamic limit lut 1, it sufhces to set fc = 
and expand the integrals in powers of lot. The subtleties 
and the typical calculation details are give in Ref. [THI l. 
to which we will refer in what follows. 

The power attenuation can be obtained from 



where 



where 



a{uj) = -2-/{uj)/vs, 



7H = i^o(fc) im [n^H], 



(34) 



(35) 



and Vg is the sound velocity. 

The renormalization of the phonon frequency a;(fc) is 
obtained from the real part of phonon self-energy using 



AiE) 



\^p\')f.s. 



F.S. 



(40) 



with Ap the superconducting energy gap, E = 
y^^ + |App, and f{E) the Fermi distribution function. 
For longitudinal phonons propagating in s- and d-wave 
superconductor, T^p does not depend on ip. For trans- 
verse phonons in s-wave superconductor the attenuation 
does not depend on the polarization direction because 
the energy gap Ap is isotropic in the plane, and because 
of Eq. Of course, the dependence of the attenua- 

tion on temperature can be different for different order 
parameter symmetry. |^ ITol . ITU 

Let us consider the fourth case, namely, transverse 
phonons in d-wave superconductor, at the energy gap 
given by Ap = Aq cos{2ipp). Let phonons propagate per- 
pendicular to the conduction layers (along the z axis. Fig. 
I2J), with the angle ip between the phonon polarization 
vector e and the x axis in plane. The electron-phonon 



Lu{k) = wo(fc)VTTRelTR. 



(36) 



For both longitudinal and transverse sound polariza- 
tions, one finds: 



1 

2, 



(37) 



where index e denotes the phonon polarization (T or L). 
One can see that for transverse polarization, the phonon 
self-energy does not depend on the orientation of the vec- 
tor e in the plane. This is true for both tight-binding 
form of the spectrum Eq. (|24|1 and low-density limit Eq. 
(I16|l . In the latter case, one can analytically obtain: 



\ip ~ (pp)d(pp 



(38) 



that does not depend on ip. 



B. Sound attenuation at T < Tc 

Before analyzing the fluctuation effects, let us consider 
the effect of the in-plane phonon polarization and order 
parameter symmetry on sound attenuation at tempera- 
ture quite below Tc. It is possible to show that the 
sound attenuation in the Born limit is given by: 



a(T)_ 



dE 



E '■ 



(39) 




FIG. 2: Sketch showing the d-wave superconducting energy 
gap. The ultrasound is propagating along k direction with 
the polarization e having an angle ip with the x axis. 



vertex is then given by Eq. H27|) . As was pointed out in 
Sec. IIV Al in the normal state the sound attenuation and 
sound velocity do not depend on the orientation of the 
polarization in plane. Combining the expressions Eqs. 
(|27|l . H39|l . and H4(J|) . one can see that the sound attenua- 
tion does not depend on the polarization direction below 
Tc as well. This can be demonstrated analytically in low- 
density limit using Eq. (|38|l . and numerically in more 
general case of a square quasi-2D lattice. This should be 
contrasted with the results of Refs. and , where 
the attenuation of transverse phonons with both k and 
e lying in the plane does depend on the propagation di- 
rection. 
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V. FLUCTUATION CORRECTIONS, S-WAVE 
SUPERCONDUCTOR 

A. Generalities 

The Feynman diagrams that give the main contribu- 
tion to the renormalization of the electron-phonon loop 
Eq. Ip^ by superconducting fluctuations are presented 
in Fig. |3| Here, each wavy hue corresponds to the fluc- 

3 






FIG. 3; Feynman diagrams that give the leading-order cor- 
rections from the superconducting s-wave fluctuations to the 
sound attenuation as T ^ Tc. Diagram 1 is of the Aslamazov- 
Larkin ( AL) type, diagram 2 is of Maki- Thompson (MT) type, 
and the diagrams 3-6 are of the density of states (DOS) type. 
Sohd lines are the normal state Green functions, wavy lines 
the fluctuation propagators, shaded semicircles the impurity 
ladder averaging, dashed lines with cross the single impurity 
scattering, and open circles the renormalized electron-phonon 
vertices. The corresponding diagrams with a d-wave order pa- 
rameter symmetry do not contain impurity semicircles (See 
Ref. [23). 

tuation propagator (Cooper ladder) L{q, flk) which, as 
T~>Tc, has the form 

L^''\q,nky^ = -uo \e + T]q^^+rsm'^{q,c/2) + nki 

(41) 

where e « (T — Tc)/Tc, flk is the bosonic Matsubara 
frequency, and the coefficient 77 has the meaning of the 
square of the effective coherence length ^ in the isotropic 
2D case for s-wave pairing: 



r^^e (Tt) = - 



(tvf) 



(42) 



Attx J 



4ttx 



4' 



where ipiz) is the digamma function. Then, r = 
IGt^rj/vp ^ 1 is the anisotropy parameter |2g. At 
T = Tc , the anisotropy parameter can be written as 
r = 4^^(0)/c^ where Ci(O) is the Cooper pair size in 
the perpendicular (z) direction. Finally. Ts = 7r/(8T) is 
the Ginzburg-Landau time. 

Returning to Fig.|3 the shaded semicircles correspond 
to vertex corrections from impurity averaging and arc 
given by: 



A(q,£i,e2) 



|£l - £2] 



El - £2! 



TQi~£l£2) 



(43) 



We neglect the contribution from the diagrams that con- 
tain the impurity ladder in the particle-particle channel 
1151 1 because they arc less singular in e (they arc not 
shown on Fig. |3J). 

The open circles at the extreme left and right-hand 
sides of the diagrams in Fig. El represent the electron- 
phonon vertices Tep{p), that contain the dependence on 
electron momentum and that will be different for longi- 
tudinal and transverse phonon polarizations. Note that 
there is no impurity averaging of the vertices in fluctua- 
tion diagrams just like in the normal state polarization 
operator in the hydrodynamic limit. 

As usual, the diagrams in Fig. |3| correspond to three 
different manifestations of superconducting fluctuations. 
[T^ a) Some of the electrons behave like Cooper pairs for 
a time given by the Ginzburg-Landau time. This is the 
famous Aslamazov-Larkin (AL) contribution Fig. Oil), 
b) The single-particle excitations are Andreev reflected 
off the superconducting fluctuations, as described by the 
so-called Maki-Thompson (MT) term Fig. 012). c) The 
effective number of normal carriers is reduced because 
some of the electrons exist as transient Cooper pairs. 
This is the so-called density of states (DOS) contribution 
in Fig. 121(3-6). Additionally, the analytical expression for 
the MT term can be separated into the anomalous (aMT) 
and regular (rMT) parts. The difference between them 
comes from the fact that there is additional diffusion pole 
in the integral for the anomalous MT part, that in most 
cases enhances the aMT fluctuation corrections. 

The evaluation of the integrals corresponding to the di- 
agrams in Fig. 13 includes taking the limit 17^ = in the 
Green's functions and also in the fluctuation propagator 
of the DOS and MT diagrams, and analytically continu- 
ing the diagrams to real phonon frequencies —i- —iuj. 
In the hydrodynamic limit, where the electron mean free 
path i is much smaller than the wavelength of sound and 
where the electron collision rate is much larger than 
the sound frequency, we expand the result in powers of 
LUT and we set the phonon momentum k ^ because 
the expansion in powers of k£ gives corrections negligible 
to leading order of cut. The terms with odd powers of 
Lu contribute to the sound attenuation while the terms 
proportional to even powers of uj contribute to the sound 
velocity. In the most general form the corrections can be 
written as: 

Aa('^-^)(T,c.) = 5|^«(f.^.^)(Tr)/(^-)(e,r,7,), (44) 
epVs 



(45) 



£f 



'(rr)/f-)(e,r,70). 



Here /3 denotes the particular channel (DOS, rMT, aMT 
or AL), superscript e the phonon polarization (L or T). 
The superscript s stands for the s-wave symmetry of 
the order parameter and in following Sections on d-wave 
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symmetry it can be also d. The symbol (e, r, 7^) 

denotes the function of temperature which usually con- 
tains the main singularity and comes from the integration 
over the momentum q in the fluctuation propagator Eq. 
(|41|l . and which also depends on the material properties. 
The 70 = 2ri / (vjpTTif,) in the temperature functions is 
the cutoff parameter that appears in the anomalous MT 
integrals and that depends on the phase-breaking time 
T0. [iJI When the subscripts a or 1; on f^^'^\e, r, 7^) arc 
omitted it is because the function is identical for the at- 
tenuation and sound velocity cases. 



B. Longitudinal phonons 

For s-wave longitudinal case, the analytical expressions 
for the temperature functions f^^'^^ read: 



/('^*^^-)(6,r,7,) 
1 



"f<p 



\/l4>il<t> + r) + y/e{€ + r) 



r 



2(e + r) 



2e , 
1 + — ( -1 

r 



(46) 

(47) 
(48) 



(49) 
(50) 



The coefficients Ki^'^'''\TT) in Eq. (gH) and Ki^'^''\TT) 
in Eq. (|45|l come from the integration of Green functions 
and impurity blocks. They have weaker temperature de- 
pendence as e than /(z^'^' (e, r, 7^), and they basically 
show how the impurity concentration affects the result. 
In Eqs. H49|l and (|5()|) we introduced indices a and v to 
distinguish between the temperature corrections to the 
sound attenuation and the sound velocity that are differ- 
ent in the AL channel. 

The analytical form and an extended discussion of the 
asymptotics of the /''^'''^(e, r, 7^) functions and the k 
coefficients can be found in Ref. In short, the 

temperature functions experience sharp enhancement as 
e ^ 0, although they remain finite at e = 0. Note that 
j(rA/T,.)(g o,r) = /i^^'''(e = 0,r) = 1, while the 
value of functions /(^°^''*)(e, r), /("^^^'^' (e> 70) and 
fj^^'''\e,r) can be much larger than 1 at high mate- 
rial anisotropy. In addition, long phase breaking times 
increase Z'"*^"^'*-' (^i '^i T^)- Thus, with the parameters 
corresponding to the real materials, the rMT diagram is 
always negligible. In contrast to the conductivity fluctu- 
ations, the contributions of the DOS diagram is enhanced 
and becomes comparable to that of the AL and aMT di- 
agrams, although its sign is opposite to that of AL. The 
superconducting fluctuation corrections to the sound at- 



tenuation, in a realistic range of parameters, are given 
by the sum of the DOS, anomalous MT and AL terms, 
which decrease the normal state attenuation. In contrast, 
to leading order in cu, the corrections to the sound veloc- 
ity are given by DOS diagrams, because the expansion 
for the anomalous MT diagram begins at order uj while 
the AL diagram is small to this order. 

One can also see that the temperature dependence of 
the microscopic terms Eqs. (|46|) . (|47f) . and (|50|l is in 
agreement with the thermodynamic results Eqs. H19(l - 
(|21|l that depend on the modulation of the hopping inte- 
gral. 

Summarizing, the signs of the principal terms for s- 
wave longitudinal case is as follows: 

and all other terms can be neglected. 



C. Transverse phonons 

In the s-wave transverse case, the structure of the DOS 
and MT terms is as follows: 



(2^ 



(52) 

L(«)(g,l]fe)if(^'^'^)(g,17,,c..), 



where K^^''^'"'' {q,Qf;,uj^) is the "bubble" that contains 
electronic Green functions Eq. (|33|) . impurity vertices 
Eq. H43|l . and electron-phonon vertices Eq. H27() . There 
is no significant modification of the integrals in compar- 
ison with the longitudinal polarization case. Indeed, the 
angular dependence of the electron-phonon vertex part 



in the DOS diagram is averaged out at the 

Fermi-surface. 

In the Maki-Thompson diagram, the term 
Fep (tp, p)rep ((/3, q — p) after averaging over the 
direction of p provides an additional term cos{qxa) in 
the integrals over q: 



d'^q cos{qxa) cos{qzc) 
(27r)3 e + r/qf, + r sin^(q^c/2) ' 



for the regular part, and 



dh 



cos 



i{qxa) cos(g^c) 



- r sm 



'■{q.c/2) 



(53) 



(54) 



for the anomalous part of the diagram. (We set ^3 = 
from now on, assuming the polarization is along the a;- 
axis. The generalization of our results to arbitrary po- 
larization is straightforward and does not change our 
main conclusions). At app <C 1, one can neglect the 
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(72:-depcndcncc in cos(qxa) in the above integrals and im- 
mediately obtain the longitudinal MT terms Eqs. H47|) 
and (|49|1 . with the properly defined transverse coupling 
constant gx- 

The Aslamazov-Larkin diagram is given by: 

xL^-^Hq,nk)L'^-'Hl,^k+io,), (55) 

where B{q,flk,i-u^) is the triangular block in the dia- 
gram (1) on Fig. 13 As usual, at g — s- 0, which gives 
the main contribution from fluctuation propagators, the 
Green function G{q — p) in the B-block should be ex- 
panded in powers of q: 

G{q -p) = G{-p) + G'{-p)Aaq,p) (56) 

+ G'{-p)Aaq,p)' + ... 

where A^{q,p) = ^{q — p) — ^{p)- One can see that with 
the quasiparticle spectrum Eq. H16|l one obtains 



Ae(g,p) 



2m 



(57) 



2tj_{cos[((7^ ~Pz)c] - cos{pzc)}. 



The first term of the expansion Eq. H56() that does 
not vanish upon angular averaging is G(~p)^{AS,{q,p))'^. 
Thus the B-block contains a product of five Green func- 
tions; one can check that this does not lead to any sin- 
gularity in the total AL diagram contribution. 

We must also consider the term in the quasiparticle 
energy spectrum that depends on hopping tx t± along 
the diagonal lattice bonds R = ±a±c (See Fig. 0). This 
term can be chosen as: 



such that 



^^iq,p)x 



i{p)x = tx sinip^a) sin(p^c). 



qxjqx - '2px) 
2m 



(58) 



(59) 



+tx{fiin[iqx - Px)a\ sin[((7^ - p,)c] - sinip^a) sin(p^c)}. 

With this modification, the first non-vanishing term 
of the expansion Eq. (|r)6|l appears to first order in 
A£^{q,p)x- The integration in the B-block now is reduced 
to J {repA^)p,s.G{p)G{p+k)G^{—p)d£,p; one can see that 
the integration over d^p gives the same result both for 
longitudinal and transverse polarizations. Indeed, there 
is no angular dependence in the integral / G{p)G{p + 

k)G^{—p)d^p. Then, the average (rip' A^a')_f.5. contains 
a factor: 



grtx — [1 - cosjaq^) cos(cg^)]. 



while its analog in the longitudinal case (rip' A^^) p.s. is: 



(60) 



Integrating over q in Eq. (|55|l with Eq. (|60|l one can set 
cos{aqx) ~ 1, thus reducing it to Eq. H61(l . From Eq. 
(26) of Ref. one can easily obtain the expression for 
the B-block: 



B{q) 



txWo[l - cos{cqz)] 



(62) 



We see that the difference between the longitudinal and 
transverse sound attenuation and velocity in AL diagram 
is only in the prefactor. However, compared with the 
DOS and MT contributions, this prefactor contains an 
additional {tx/ti_Y term because of the anisotropy pa- 
rameter r which enters the fluctuation propagator which 
is still proportional to t\. One would expect that this 
prefactor satisfles {tx/t±)'^ ^ 1, thus making the AL 
term negligible. 

In summary, the main contribution to the fluctuation 
corrections is given by the DOS diagrams (and aMT di- 
agrams in sound attenuation at sufficiently long pair- 
breaking times). The temperature function of the rMT 
diagram is less singular and the AL diagram contains an 
additional small factor, thus being negligible. The signs 
of the principal terms are as follows: 

^^{DOS,T,s) < 0^ ^aiaMT,T,s) ^ (gg) 

^^iDOS,T,s) > 0^ 



VI. FLUCTUATION CORRECTIONS, L'-WAVE 

A. Modification of the d-wave integrals 

If we assume that the virtual Gooper pairs have d- 
wave symmetry, the modification of the calculations is as 
follows. 

We have to take into account the momentum depen- 
dence of a pairing interaction that now can be written 
as: 



V{p,p') = 7]{'fp)gr]{ipp,), 



(64) 



where, for d-wave pairing, 'r]{ipp) ~ cos{apx) — cos(apy) cx 
\/2 cos{2tpp) . 

Then, in all fiuctuation diagrams of Fig. |31 each wavy 
line corresponds to the fluctuation propagator with ad- 
ditional symmetry factors: 



r]{(pp)L^'^\q,nk)f]{(pp'), 



(65) 



gLt± — [l - COs{cqz)]. 



(61) 



where L^'^\q, Vlk) has the same form as the s-wave fluc- 
tuation propagator Eq. (|41|l with modified e, ^, r, and 



Ld{q, ^l) ^ = -vq ed + Vdqf\ + Td sin^(cg^/2) - iVlTd 

(66) 

In the quasi-2D case, and assuming that the phase- 
breaking processes contribute only as the upper cutoff 
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for the Maki-Thompson diagram, the coefBcicnts in Eq. 
(|66|l read (compare with those of Ref. H^): 



e,i = hi 'ih I — 



1 



1 



'1 



1 - — V^' ( - + — ^ 

47ra; \ 2 47rx / 



, , 1 1 
^ ' 2 47rrT 



1 



(67) 




\8ttx) 



1 



1 



^"'2 + 4^) 



(68) 



f-V 

V2T/ 



x < 1, 



7C(3) ^^1. 



^ "i^'^' + 4^) 



(69) 




(70) 



X > 1. 



where a: = Tr. Note that the dependence of rj^ and r^; 
on the impurity concentration (through the parameter 
Tt) is the same. It also should be noted that the d-wave 
parameters ed, rjd-, Td, and Td coincide with the corre- 
sponding s-wave terms in the clean limit (Tr 3> 1). In 
the opposite, dirty limit, the pair-breaking effect of impu- 
rities for d-wave pairing makes these coefficient different 
from those for s-wave pairing. 

Note also that in the expression Eq. (|?)7jl Tc is the real 
superconducting transition temperature. In the case of 
d-wave pairing, it is renormalized even by scattering on 
non-magnetic impurities, and it satisfies the Abrikosov- 
Gorkov equation: 



1 



= 0. 



In the clean case, at tT w tTc ^ 1, there is not much 
difference between Tc and Tco- In the opposite, dirty 
limit, when the elastic scattering time is small, d-wave 
superconductivity is suppressed and Tc/Tco ^ as tT 
0. On the other hand, there is no such impurity effect on 
the s-wave Tc. 

Finally, there is no impurity averaging in the fluctua- 
tion propagator (no shaded semicircles on diagrams Fig. 
13). [231 This also leads to the absence of an anomalous 
part in the MT diagram 2 of Fig. |21 since there is no 
diffusion pole in the corresponding integral. 



B. Longitudinal phonons 

1. DOS 

The corrections to the sound attenuation and to the 
sound velocity have the structure of Eqs. H44|) and H45|l . 
with the temperature functions appropriate for d-wave 
pairing: 



(71) 



f-d + \Ji-d + Td 

lnr<j , ed < Td, 
In Ed , Td < ed, 



with ed and Td given by Eqs. (|67|l and l|70() . The renor- 
malization of Td Eq. (|69|l does not affect the results be- 
cause the main contribution from the fluctuation propa- 
gator is at — > 0. 

The f^(DOS.L,d) coefficients are inversely proportional 
to the renormalized rjd = in Eq. (|68|l which is dif- 
ferent from its s-wave analog. Another source of differ- 
ence between s-wave and d-wave results is the absence of 
impurity vertices A in d-wave diagrams. With all these 
modifications, one finally obtains the dependence of the 
i^(DOS,L,d) coefficients on sample cleanness through the 
parameter Tt « TcT. 

For sound attenuation^ the coefficient is given by: 



ADOS,L,d) 



1 







a- 


4nx} 



2 Attx J 
4(27r-3)x, a;<l, 

/I ^3^2 

X > 1, 



+ 8n{2n ~ l)x4'" - + 



7C(3) ' 

while for sound velocity it reads 
mos,L,d) ^ 



(72) 

Anx J 



(73) 



(74) 



The temperature function f^^^^''^\ed,rd) Eq. lf7T|l is 
similar to its s-wave analog Eq. (|46|) . although now it 
also depends on the impurity concentration. Qualita- 
tively, it is monotonically increasing as T —t Tc just like 
/(^0^-^)(e,r). 



2. rMT 

Modifications similar to those in the DOS integrals also 
occur in the MT diagram. There is no aMT contribution 
because there is no diffusion pole originating from the im- 
purity vertices. For rMT term, the temperature function 
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1 , 

Id_ 

4ed 



rd 

< fid- 



(75) 



Similarly to the s-wave case, the function in Eq. (|75|l has 
a limit /^''^^^'''Herf = 0,^^) = 1. This fmiction is very 
little affected by Tr. The ni^^TXA) coefficients are finite 
(they do not contain a singularity as functions of x). In 
summary, the rMT contributions to the attenuation and 
sound velocity is negligible in comparison with those from 
the DOS diagrams. 



AL 



The triangle of Green functions in diagram 1 in Fig. El 
does not contain impurity semicircles and reads: 

xG{p, e„)G{p, e„ + uj„)G{q - p, -£„). 

One can see that just like in the "longitudinal s-wave" 
case, it suffices to expand the Green function G{q — 
p,—£n) to first order in hopping ti_. The evaluation of 
Eq. (|75|l with ri{ipp) = ^/2cos{2tpp) gives the same result 
as in the case of s-wave pairing with r]{ipp) = 1. Then, 
the temperature functions of the AL diagram can be ob- 
tained from those for the s-wave pairing Eqs. H49|) and 
(|50|l with the substitution of e^, rjs, t^, and by their 
d- wave analogs Eqs. (|^ - (|7nil . 

In the clean limit we recover for K^f^^'^''''^ and kI"^^'^'''^ 
their longitudinal s-wave, Tr-independent analogs [T5l |. 
On the other hand, i^i^^-'^'d-) coefficients remain finite in 
the dirty limit. Moreover, in highly anisotropic materials 
at r^i <^ Erf, the AL contribution to the sound attenua- 
tion is clearly suppressed, and becomes important only 
at temperatures extremely close to Tc, just like the corre- 
sponding s-wave terms. The AL correction to the sound 
velocity is negligible just like the rMT one. 

4- Summary 

In summary, in the d-wave fluctuation corrections to 
longitudinal sound the aMT term is absent, and rMT 
and AL sound velocity corrections can be neglected in 
the range of parameters that is reasonable. The signs of 
the principal terms arc as follows: 

^^iDOS,T,s) > 0^ 



The magnitude of the corrections is finite though large 
in the clean limit, and it is suppressed to zero in the 
very dirty limit. Again, like in the s-wave case, the tem- 
perature dependence of the microscopic expressions for 
d-wave order parameter symmetry can be obtained in 
the phcnomenological approach with the substitution of 
proper Ed, ??d and in the fiuctuation term Eq. (|17() and 
taking into account the strain dependence of t± . 



C. Transverse phonons 

One can easily combine the results of Sec. IV Bl and 
Sec. rvTBl 

The DOS contribution is given by Eq. (|52() with the 
fiuctuation propagator modified according to Eq. H65(l . 
It is clear that after such a substitution one should obtain 
the results for the DOS longitudinal d-wave term because 



the angular dependence in 



ep 



is averaged out. 



Similar arguments are valid for the rMT contribution. 
Indeed, the longitudinal term contains the combination 



r(^)(p)ri^)(q- 



p) 



(78) 



9l cos(p^c) cos[(q^ - p^)c]. 



The corresponding combination in the transverse term 
at, for example, ip ~ can be rewritten to first non- 
vanishing order in the low-density limit app <^ 1: 

rip(^ = o,p)r(p(^==o,q-p) 



(79) 

= sin(p^a) sin(pj,c) sm[{qa: - Px)a] sm[{q^ - p^)c] 
« cosiq^a) cos^((pp) sm{p^c) sm[{p^ - q^)c], 

where (app)^ is absorbed in gx, and the term linear in 
cos(iy9p) drops out after angular integration. 
Since 



7r/c 



cos(p^c) coa[{qz - Pz)c]dp^ 



(80) 



— tt/c 



/ sin(p^c) sin[(q2 - Pz)c]dp^ = - cos(g^c). 



one can see that after integrating over pz , the transverse 
term should contain an extra cos{qxa) cos^ {(pp) in com- 
parison with the longitudinal one. Then, it suffices to 
set cos(ga;a) « 1 at qxO ^ 1, and cos^{ipp) will provide 
an extra 1/2 after the angular integration. We conclude 
that with the proper choice of the coupling constant, one 
can use the longitudinal d-wave results. 

For the same reason as in the longitudinal d-wavc case, 
there is no aMT contribution here. 

Finally, as described in Sec. IV Bl in order to obtain 
non-zero AL contribution, one should expand the Green 
function G{q—p) either to second order in small A^{q,p)± 
as in Eq. H57|l (this results in non-singular AL contri- 
bution), or to linear order in A£^{q,p)x (this adds an 
additional small factor {tx/t±f' in front of the whole 
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diagram). The details of the integration over q are the 
same as in the transverse s-wave case. We conclude that 
the AL term can be neglected too. 
The leading terms then read: 

^^{DOS,T,s) 0^ ^y{DOS.T,s) > ^ (g^) 

VII. DISCUSSION 

The superconducting fluctuations can provide a renor- 
malization of the normal state sound attenuation and 
sound velocity that depends on the various phonon po- 
larizations and order parameter symmetries. The mag- 
nitude of the effect, as can be seen from Eqs. H37(l and 
(|44|l . is of the order of Tc/ep or less, depending on the 
material cleanness. For layered organics, the Fermi sur- 
face parameters of Ref. |23 lead to Tc/ep ~ 10~^. The 
temperature functions /'•''''^^(e, r, 7^) increase this ratio 
at T ^ Tc, thus making the fluctuation corrections ex- 
perimentally measurable. 

The actual sign of the fluctuation corrections is deter- 
mined by the microscopic coefficients k^'^''^''''. We gather 
the information on signs and relative magnitudes of the 
microscopic fluctuation terms in Table I. In short, the 
leading contribution is always in the DOS channel, and 
there are comparable corrections from the aMT diagram 
only for s-wave attenuation but for both polarizations, 
and from the AL diagram only for longitudinal sound at- 
tenuation but for both symmetries. The rMT (always) 
and the AL (longitudinal sound velocity) temperature 
functions have a limit unity at T = Tc but could com- 
pete with the leading terms only in the special cases of 
low anisotropy and short phase breaking time where the 
leading terms are suppressed. The aMT terms in the 
s-wave sound velocity and the AL terms for transverse 
phonon polarization contain additional small factors and 
can thus be neglected. Finally, there is no aMT term 
for d-wave symmetry of the order parameter because the 
diffusion pole disappears in this case. 

The temperature dependence in the sound velocity ob- 
tained from the phenomenological Ginzburg-Landau free 
energy is in agreement with the microscopic calculations 
for the longitudinal phonons in both the s-wave and d- 
wave cases. In the latter case, it suffices to replace e, ry 
and r by their d-wave equivalents e^, 77^ and defined in 
Sec. IVII For example, one can compare the temperature 
dependence of the phenomenological result Eq. (|19|l with 
that of the DOS contribution in Eqs. and (TI])- Also, 
the temperature dependence of the phenomenological Eq. 
(|20|l as well as the sign of the correction is the same as 
that of the rMT terms Eqs. 1071) and (|7SJ). Similarly, the 
AL contributions (sign and temperature function) Eqs. 
(|50|l and Sec. IVI B 31 arc explained by the phenomeno- 
logical result Eq. H21|) . The Ginzburg-Landau approach 
does not treat the aMT contribution correctly because it 
does not include the dynamics involved in the diffusion. 

In principle, one can obtain the fluctuation correc- 



tions to the sound velocity of transverse phonons from 
the Ginzburg-Landau approach as well. The modiflca- 
tion will include taking into account the next-to-nearest 
neighbor hopping in the quasiparticle energy spectrum 
Eq. (|16|l . extra hopping terms in the fluctuation free en- 
ergy Eq. Hll|l . and choosing the appropriate combination 
c in Eq. Q that will modify the elastic terms Eq. (j^J). 

It should be noted that the relation cx — be- 
tween fluctuation corrections to the sound velocity and 
specific heat that is valid in the bulk continuous model, 
is not satisfied in the tight-binding model with hop- 
ping terms in the electron quasiparticle energy spectrum. 
Thus, the experimental information from the specific 
heat measurements (first reference in 0) is not enough 
to predict the behavior of the sound velocity fluctuation 
corrections in organic compounds. 

We have obtained these results in the quasi-2D square 
lattice model for a corrugated cylindrical Fermi-surface. 
As a result, we see, for example, that there is no ef- 
fect on the attenuation of the transverse sound as the 
direction of polarization is changed in the plane, both 
in the normal and in the superconducting states. In or- 
der to make a better fit of the experimental data, the 
generalization of our model to the particular crystal lat- 
tice (that is not tetragonal) and to the appropriate band 
structure is important. For example, it might be nec- 
essary to incorporate the dimcr structure and triangular 
in-plane symmetry of organic materials. |2^|33,|^ That 
would better refiect the two-band energy spectrum of the 
K-(ET)2X compounds. All these problems should be the 
subject of future research. 

In conclusion, we have found how the sound atten- 
uation and sound velocity are renormalized by super- 
conducting fiuctuations at temperatures close to Tc in 
layered superconductors. For various polarizations of 
phonons propagating perpendicular to the conduction 
plane, we considered the cases of s- and c?-wave symme- 
try of the order parameter. In the hydrodynamic limit 
LOT <^ 1 and fci? <C 1, we found the contributions from 
all the fluctuation diagrams (namely, the AL. MT and 
DOS) and provided the theoretical background for the 
analysis of the experimental information in organic su- 
perconductors. The complete analytical expressions for 
the microscopic fluctuation corrections can be found else- 
where • 
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TABLE I: Summary of the fluctuation contributions from different diagrams for various order parameter symmetries and 
phonon polarizations. The symbol < | — 1| means that the sign of the correction is negative and that the absolute value of the 
temperature function is less than unity with limit unity a.t T = Tc- 
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